Phonon-like and single particle dynamics in liquid lithium 
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The dynamic structure factor, S(Q,E), of liquid lithium (T—475 K) has been determined by 
inelastic x-ray scattering (IXS) in the momentum transfer region (Q— 1.4-^110 nm _1 ). These 
data allow to observe how, in a simple liquid, a phonon-like collective mode evolves towards the 
single particle dynamics. As a function of Q, one finds: i)at low Q's, a sound mode with a positive 
dispersion of the sound velocity, it) at intermediate Q's, excitations whose energy oscillates similarly 
to phonons in the crystal Brillouin zones, and Hi) at high Q's, the S(Q,E) approaches a Gaussian 
shape, indicating that the single particle dynamics has been reached. 



PACS numbers: 61.25.Mv, 61.10.Eq, 61.20.Ne 

The spectrum of density fluctuations, S(Q,E), of a 
liquid shows a very rich phenomenology with features 
strongly dependent on the considered momentum, Q, 
and energy, E, regions The shape of the S(Q,E) 

is well established in the small and high Q limits. In 
the hydrodynamics limit, i. e. at Q/Qm << 1 (Qm is 
the position of the first maximum of the static struc- 
ture factor S(Q)), the spectrum shows three peaks - the 
Brillouin triplet. They are respectively the Stokes and 
AntiStokes propagating compressional modes -dispersing 
linearly with the adiabatic sound velocity v s , and the 
thermal diffusion mode centered at zero energy transfer. 
In the opposite limit, i. e. at Q/Qm >> 1, one reaches 
the so called impulse approximation, where the excited 
particles acquires a kinetic energy much higher than any 
inter-particle potential energy. Therefore the target par- 
ticle behaves as a free particle, and the S(Q,E) line- 
shape reflects the particles' initial state momentum dis- 
tribution. Considering the Boltzmann distribution, the 
S(Q, E) reduces to a Gaussian centered at the recoil en- 
ergy h 2 Q 2 /2M, and with variance H 2 K B TQ 2 /M. Here, 
M is the particel mass. The evolution between these two 
limit cases is affected by a realm of dynamical processes 
as the interaction of sound waves with other degrees of 
freedom (translational diffusion, rotations and internal 
modes for molecules), and the interaction between dif- 
ferent collective modes, responsible, for example, for the 
slowing down of the diffusional dynamics in supercooled 
liquids. Moreover, when Q is comparable to Qm, impor- 
tant modifications of the Brillouin triplet occur also as 
a consequence of structural effects. In fact, the sound 
waves in the liquid cannot be longer described as density 
fluctuation of a continuum medium, and the local struc- 
ture with its intrinsic disorder becomes relevant. Com- 
pared to the crystalline case, the absence of long range 
order will prevent to a certain extent, in the liquid, the 
replication of the sound dispersion relation in high order 
Brillouin Zones (BZ). 

This rich phenomenology has motivated, since long 



time, the experimental study of the dynamics of liquid 
systems. Using ultrasound absorption methods and Bril- 
louin light scattering techniques, the sound waves, and 
their interactions with the relaxation processes active 
in the liquid, have been studied in great detail up to 
Q/Qm ~ 10~ 3 . Similarly, Inelastic Neutron Scattering 
(INS) has been used at very large Q transfers to deter- 
mine the S(Q, E) in the impulse approximation, and in Q 
regions close to Qm to determine the S(Q,E) lineshape 
when the two limit description are expected to fail ||Q . 
The kinematic limitations on the accessible Q — E region 
of existing INS instruments || did not allow to study 
with continuity and under comparable experimental con- 
ditions the liquid dynamics in the Q range spanning from 
the collective to the single particle behaviors. Moreover 
the neutron scattering cross-section accounts for two dif- 
ferent contributions: besides the coherent cross section 
-probing the collective dynamics- there is an incoherent 
contribution which, at each wavevector, reflects the single 
particle motion and therefore hides the crossover between 
the previously quoted regimes. The devolepment of In- 
elastic X-ray Scattering (IXS) § has allowed, recently, to 
sensibly extend the accesible Q — E region in disordered 
materials and to avoid the incoherent contribution. 

In this work we report the determination by IXS of the 
dynamic structure factor of liquid lithium in the 1.4-^110 
rnir 1 Q range, corresponding to Q/Qm ~ 5 ■ 10~ 2 — 
5. Lithium has been chosen because, among the simple 
monoatomic liquids, is the one that is best suited to IXS. 
Indeed, its low mass gives recoil energies observable in the 
considered Q-range, and its low atomic number and large 
sound velocity («5000 m/s) give optimal signal with the 
available energy resolution, compensating for the large 
form factor decrease at high Q values jjj. The S(Q,E) 
spectra, reported to their absolute scale exploiting the 
zcroth and first moments sum rules, show the transition 
from a triplet to a Gaussian. The maxima, Q(Q), of 
the longitudinal current spectra (E 2 j Q 2 S (Q , E)) show 
an almost linear dispersion relation at low Q - typical 
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of a sound wave - and a completely different dependence 
in the high Q limit, where it approaches the parabolic 
dispersion of the free particle. These two regions are 
joined by oscillations of f2(Q), which are in phase with the 
structural correlations, as observed in the S(Q). Finally, 
in the low Q region these data confirm the existence of 
positive dispersion of the sound velocity. 

The experiment has been carried out at the very high 
energy resolution IXS beamline ID16, at the European 
Synchrotron Radiation Facility. The incident x-ray beam 
is obtained by a back-scattering monochromator operat- 
ing at the Si(h h h) (h=7, 9, 11) reflections g. The 
scattered photons are collected by spherical silicon crys- 
tal analyzers, operating at the same Si(h h h) reflection 
The total energy resolution — obtained from the 
measurement of S(Qm, E) in a Plexiglas sample which is 
dominated by elastic scattering — is 8.5 meV full-width- 
half- maximum for h=7, 3 meV for h=9 and 1.5 meV for 
h=ll. The momentum transfer, Q = 2khsin(9 s /2), with 
kh the wavevector of the incident photon and 9 S the scat- 
tering angle, is selected either between 1.4 and 25 nm" 1 
by rotating a 7 m long analyser arm in the horizontal 
scattering plane (data taken at h=9 and 11), or between 
24 and 110 nm" 1 by rotating a 3 m long analyser arm in 
the vertical scattering plane (h=7). The total Q resolu- 
tion has been set to 0.4 nm -1 at h=9 and 11, and 1 nm -1 
at h=7. On the horizontal arm, five independent analyz- 
ers were used to collect simultaneously five different Q 
values, determined by the constant angular offset of 1.5° 
between neighbour analyzers. The vertical arm houses 
only one analyzer. Energy scans are done by varying 
the temperature of the monochromator with respect to 
that of the analyzer crystals. The absolute energy cal- 
ibration between succesive scans is better than 1 meV. 
Each scan took about 180 min, and each (Q-point spec- 
trum has been obtained from the average of 2 to 8 scans 
depending on h and the Q-transfer. The data have been 
normalized to the intensity of the incident beam. The liq- 
uid lithium uncapped container is made out of austcnitic 
stainless steel with a resistence heater, used to kept the 
liquid at T = 475 K. The 20 mm long sample, kept to- 
gether by surface tension, was maintained in a 10 -6 bar 
vacuum. The lithium has been loaded in an argon glove 
box. In the Q—E region of interest, empty vacuum cham- 
ber measurements gave either the flat electronic detector 
background of 0.6 counts/min or, at 9 < Q < 13 ran -1 , a 
small elastic line due to scattering from the chamber kap- 
ton windows (each 50 mm thick) which was subtracted 
from the data. 

The IXS spectra of liquid lithium are reported in Fig. 1 
at the indicated Q transfer values. The low Q data show 
the Brillouin triplet structure with the energy of the in- 
elastic peaks increasing with Q up to a Q value of 12 
nmT 1 . This value corresponds to Qm/2, as deduced from 
the S(Q) reported in Fig. 2. One can interpret, there- 
fore, the dispersion up to Qm/2 as that of the longitudi- 



nal acoustic branch in the pseudo-first BZ. Furthermore, 
similarly to what it is found in the second BZ of a crystal, 
we observe that, also in liquid lithium, the energy of the 
acoustic modes decreases with increasing Q from Qm/2 
to Qm- Increasing Q above Qm, i- e. in the "third" or 
higher BZs, one finds that the spectrum gets incresingly 
broader and distinct peaks are no longer observable. At 
the highest Q- values one finds that the S(Q,E) becomes 
a symmetric peak centered at energies larger than E = 0. 
Beside the observation of dispersion in a first and second 
pseudo-BZs, it is also important to note that the broad- 
ening of the excitations monotonically increase with Q to 
the extent that towards the end of the second pseudo-BZ 
the inelastic features are no longer showing a well defined 
peak. 

The previous qualitative description can be made sub- 
stantially more quantitative considering that the dy- 
namic structure factor S(Q,E) can be derived from the 
measured intensity, I(Q,E), using the zeroth and first 
moment sum rules for S(Q,E): 

mf 0) = J S(Q,E)dE = S(Q), (1) 

mf 1} = J ES(Q,E)dE = h 2 Q 2 /2M. (2) 

Considering that I(Q,E) = a(Q) J dE'S{Q,E')R{E - 
E'), where R(E) is the experimental resolution function 
and a(Q) is a factor taking into account the scattering 
geometries - efficiencies and the lithium atomic form fac- 
tor, the first moments of the experimental data, iri>f \ and 

m^-p and those of the resolution function, and m^y 

are related to and mfu by: 

m( 0) = a(Q)mf 0) mf 0) , (3) 

= a (Q)( m (Q) m (l) + m fl) m (Q))- (4) 
From the previous equation one derives that 

S{Q) = W (m (i)Ho) - mfa/mfa)- 1 . (5) 

This procedure has been utilised to put the S(Q,E) on 
an absolute scale using the experimentally determined 
I(Q,E) and R(E). The reliablitiy of this procedure is 
shown in Fig. 2 where we obtain an excellent agreement 
between the S(Q) values obtained by Eq. 4 and those 
derived by Molecular Dynamics (MD) simulation pj. 

The possibility to express S(Q,E) in absolute units 
allows to compare the IXS data with the Gaussian line- 
shape expected for the S(Q,E) when the single particle 
limit is reached. This comparison is reported in Fig. lb, 
where each solid line represents the sum of two Gaus- 
sians: 

G(Q,E) = -}=Yl -^-(E-E^^af (6) 
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where Cq = 0.08 and Cj = 0.92 are the natural abun- 
dances of the 6 Li and 7 Li isotopes, Ei = h 2 Q 2 /2Mi their 
recoil energies and af = h 2 K B TQ 2 /M l . The progres- 
sively better agreement between the data and G(Q, E) 
with increasing Q testifies the evolution towards the sin- 
gle particle behaviour, which seems to be reasonably well 
reached at the highest investigated Q values and at the 
considered temperature. 

The dispersion relation of the energy fl(Q) of the in- 
elastic signal observed in Fig. 1 is obtained by determin- 
ing the maximum of the current spectra. This allows us 
to estimate ft(Q) independently from any specific model 
for the S(Q,E). Examples of current spectra, obtained 
from the spectra in Fig.l, are shown in Fig. 3. The values 
of 0(<3) have been determined by a parabolic fit of the 
maximum region in the Stokes side. The values of £l(Q) 
are reported in Fig. 4. In Fig. 4a are also reported the 
linear dispersion of the adiabatic sound velocity (dotted 
line) and the parabolic dispersion expected for the sin- 
gle particle dynamics. We observe that f2(Q) is close to 
the adiabatic sound mode at the lowest Q values, and it 
shows a positive dispersion before reaching the maximum 
of the first pseudo-BZ. This is emphasized in Fig. 4b, 
where the low Q region is expandend, and the inset re- 
ports directly the sound velocity v(Q) = fl(Q)/Q. This 
behaviour confirms previous MD Jn|-[l3|] and experimen- 
tal |l4| data in a similar Q region on lithium and other 
alkali liquid metals. With increasing Q values, the points 
in Fig. 4a show not only a second pseudo-BZ, as already 
pointed out in Fig. 1, but also a series of oscillations that 
dump out with increasing Q - here, Q(Q) approaches the 
single particle line. These oscillations are in anti-phase 
with the oscillations found in the S(Q) (see Fig. 2) and 
can, therefore, be associated to the local order in the liq- 
uid. It is of great interest to be able to follow all the 
way from the sound mode to the free particle regime the 
evolution of £l(Q). Beside the importantce of a unified 
picture in such a wide Q region these data provide the 
workbench for a quantitative theoretical analysis of the 
shape of the S(Q,E). 

In conclusion, using IXS, we have been able to measure 
the collective dynamics in a simple liquid from the collec- 
tive modes regime, dominated by sound- like excitations, 
all the way towards the impulse approximation regime, 
dominated by single particle dynamics. This textbook 
result provides a picture on how in a simple liquid the 
structural correlations induce important deviation from 
the "continuum" picture used in the hydrodynamic limit, 
and gives effects qualitatively similar to those as the BZs 
in crystals. 

We acknowledge valuable help of H. Mueller from the 
Chemestry Laboratory at ESRF for his technical assis- 
tance during the sample manipulation. 
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FIGURE CAPTIONS 

FIG. 1 - Dynamic Structure Factor (Density Correlation Function) 
of Lithium at T = 475K measured by IXS and normalized 
according to the procedure discussed in the text. The spec- 
tra on the left column are taken with a resolution AE = 3 
meV using Si (9 9 9) reflection; Those on the rigth column 
with AE = 8.5 meV using the Si (7 7 7) reflection. In the 
latter column the comparison with a gaussian lineshape (full 
line) expected in the free particle limit is also reported. 

FIG. 2 - Test of normalization accuracy: the S(Q) values as ex- 
tracted from IXS normalization (open symbols) are reported 
toghether with computer simulation data (full line) JlC| . The 
inset shows an enlargement of the small Q region and the 
dotted line is the Q — » limit expected from the isothermal 
sound velocity. 

FIG. 3 - Current Correlation Function (J(Q, E)) obtained from 
the S(Q,E) spectra reported in Fig.l as J(Q,E) = 
E 2 /Q 2 S(Q,E). 

FIG. 4 - Sound dispersion of Lithium. Fig 4a - The energy posi- 
tion of the Stokes Peak of the Current Correlation Function 
(open cirlce) are reported toghether with the dispersin ex- 
pected in two limiting cases: low Q linear dispersion with 
hidrodynamic velocity ((dotted line) and high Q (full line) 
parabolic dispersion of the impulse approximation. These 
two dispersion curves have been compute assuming the pres- 
ence of both 6 Li and 7 Li in their natural abundances (see 
text and Eq. ^ . 

Fig 4b - Detail of the 'small Q' (below the First Sharp 
Diffraction Peak) region. The deviation of the acoustic 
sound velocity from the adiabatic value (dotted line) at in- 
creasing wavevectors is observed (positive dispersion). 
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Fig. 2 - T. Scopigno et al. 
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Fig. 4a - T. Scopigno et al. 
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Fig. 4b - T. Scopigno et al. 



